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Absorption  Properties  of  Particulates  1 

in  Arctic  Pack  Ice  and  Snow:  Influence 

on  Radiative  Energy  transfer. 

Introduction 

Observations  (Grossi  et  al,  1987),  measurements  (e.g. 

Palmisano  et  al,  1987;  Perovich,  1993),  and  radiative 
transfer  models  (Arrigo,  et  al  1992;  Grenfell,  1992; 

Wiscombe  and  Warren,  1982)  have  shown  that  particulate 
material  contributes  substantially  to  the  optical  properties 
of  sea  ice  and  snow. 

From  anecdotal  observations  (Gosselin,  pers.  comm.;)  and  a 
number  of  systematic  studies  (Pfirman  et  al,  1989;  Reiminetz 
et  al,  1993),  it  is  readily  apparent  that  drifting  pack  ice 
within  the  central  Arctic  basin  contains  substantial  biogenic 
and  lithogenic  particulate  loads.  The  role  of  sea  ice 
particulates  in  radiative  energy  transfer  and/or  related 
biologic  and  thermodynamic  processes  in  these  regions  remain 
relativly  unknown. 

During  the  1989  Coordinated  Eastern  Arctic  Experiment 
(CEAREX) ,  we  investigated  the  vertical  distribution  and  light 
absorption  properties  of  particulates  in  Arctic  pack  ice. 
Concurrent  assessments  of  the  spectral  irradiance  fields 
(above  and  below  the  ice)  were  made  in  order  to  estimate  the 
effects  of  the  particulate  matter  on  the  apparent  optical 
properties  of  the  ice  cover.  From  these  analysis,  the 


realized  and  potential  effects  of  pack  ice  particulates  on  2 

radiative  energy  transfer  are  examined. 

Materials  and  Methods; 

.qtudv  Area  and  Instrumentation:  Irradiance  measurements, 

including:  surface  downwelling  spectral  irradiances 
(E(a(0  +  ,X)),  surface  upwelling  spectral  irradiances 
(Eu ( 0+ , X) ), surface  downwelling  scalar  irradiances 
(Eod(0+,PAR) )  as  photosynthetically  available  radiation  (PAR) 

(400-700nm) ,  and  under-ice  downwelling  spectral  irradiances 
(Ed(z,X))  were  made  between  April  7  to  April  22.  All 

spectral  measurements  were  made  with  a  MERlOlO 
spectroradiometer  (Biospherical  Instruments)  with  its  10  nm 
wavebands  centered  at  410,  441  455,  488,  507,  532,  550,  570, 

589,  633,  656,  and  694  nm.  Eod(0  +  ,PAR)  was  measured  with  a 
Biospherical  Instruments  model  QSR-200  Solar  Reference 
Hemispherical  Irradiance  sensor.  Both  instruments  were 
calibrated  before  and  following  CEAREX  89  experiments  at 
Biospherical  labs  in  San  Diego. 

Under-ice  measurements  were  accomplished  by  mooring  the 
spectroradiometer  on  an  90°  articulated  arm  deployed  through 
a  75  cm  diameter  hydrohole  drilled  with  a  steam  heat  drill. 

The  instrument  was  positioned,  and  leveled,  at  the  ice/water 
interface  approximately  3  m  away  from  the  hydrohole  which  was 
covered  to  avoid  stray  light  from  reaching  the  sensor.  With 
the  radiometer  still  in  position,  the  natural  snow  cover  was 
manually  removed  from  above  the  radiometer  and  the  bulk 


diffuse  spectral  attenuation  coefficients  for  both  the  sea  3 

ice  and  snow  were  calculated  according  to  calculations 
described  below. 

Spectral  albedos.  <xiX):  A  series  of  eight  incident  and 
eight  upwelling  irradiance  measurements  (spanning  a  12  hour 
period)  were  made  over  a  20  cms  of  snow  overlying  multi-year 
ice  (MY  ice) .  Surface  albedos  were  calculated  as; 

EuiO+.l) 

(A)  =  . 

Ed(0+A) 


Diffuse  Attenuation  Coefficient  for  Snow,  Kdsnw(X) ;  The 

diffuse  spectral  attenuation  coefficient  of  snow  cover 
Kdsnw(^)  (in"l)  was  estimated  according  to  an  exponential  model 

describing  the  attenuation  of  downwelling  irradiances, 


Kdsnw (^) 


In 


Ed(+snw,^) 
Ed  ( -snw,  A,) 


(1) 


H  is  the  thickness  of  the  snow  in  meters,  and  Ed(+snw,X.)  and 
Ed(-snw,A,)  are  the  under-ice  downwelling  spectral  irradiances 
(units,  |J,E  m“2  s“l)  with  and  without  snow,  respectively. 

Diffuse  attenuation  coefficient  for  sea  ice,  Khh,^o(A,)  ;  The 

bulk  diffuse  spectral  attenuation  coefficient  for  optically 
thick  sea  ice  (Kdice  (^)  )  was  approximated  from  the  following 


relationship: 


Kdice ( ^ ) 


(2) 
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1  Ed  ( z ,  A, ) 

^  *  In  - 

z  Ed(0-,X) 


(9,10)  where  Z  is  the  thickness  of  the  ice  in  meters,  Ed(Z/A,) 

is  the  downwelling  irradiance  at  the  bottom  of  the  ice  and 
Ed(0-,A,)  is  the  downwelling  irradiance  just  below  the  surface 
q£  ^2.he  ice.  Ed(0“,X)  cannot  be  measured  directly  in  sea  ice, 

but  was  estimated  as 


Ed(0-,X)  =  Ed(0+,X)  *  (1-  Ro) .  (3) 

where  R©  is  the  specular  reflection  from  the  sea  ice  surface 
which  is  estimated  according  to  Arrigo  et  al  (1992). 

Although  the  Beers-Lambert  model  is  a  derivation  of  the 
radiative  transfer  eguation  for  a  purely  absorbing  medium  and 
snow  and  ice  can  have  large  ratios  of  scattering  to 
absorption,  the  calculated  Kd(A,)  ‘s  for  both  the  snow  and  ice 

are  reasonable  approximations  of  the  asymptotic  attenuation 
coefficients  when  optical  depths  are  sufficiently  large. 

Bemnle  collections;  Two  study  sites  (HHl  and  HH2 )  were 
established  in  hummocked  multi-year  (MY)  ice  where  under-ice 
irradiances  were  measured  and  snow  samples  and  ice  cores  were 
also  obtained.  Two  additional  cores  (designated  Ll  and  L2) 
were  obtained  from  first-year  (FY)  ice  in  an  undeformed 
refrozen  lead.  Ice  cores  were  obtained  (using  a  7  cm 
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diameter  CRREL  core  barrel) ,  sectioned,  and  allowed  to  melt 
at  2°C  in  closed  opaque  containers.  Snow  samples  were 
scooped  into  clean  Nalgene  jars  and  allowed  to  melt  at  2°C. 
Visibly  apparent  discolorations  and  isolated  pockets  (ca.  5 
to  10  cm  thickness)  of  dark  particulates  were  observed  in 
several  large  75  cm  ice  cores  that  were  being  drilled  with  a 
steam  heat  drill  during  routine  camp  operations.  Samples 
from  within  these  pockets  were  also  extracted  and  allowed  to 
melt  in  closed  opaque  containers. 

Absorption  Efficiencies  of  Individual  Particles,  Oa(X,)  ; 

Ice  and  snow  meltwater  was  filtered  through  Nucleopore 
filters  (0.4  |xm  pore  size)  at  low  pressures  (5  mmHg)  . 
Particles  retained  by  the  filters  were  immediately 
transferred  to  gelatin  coated  microscope  slides  (Iturriaga 
and  Siegel,  1989)  and  kept  frozen  until  analyzed  at  USC. 
Absorption  efficiency  spectra  of  the  individual  particles 
were  determined  using  a  universal  microscope  (Carl  Zeiss) 
equipped  with  a  type  03  photometer  interfaced  to  a  tungsten- 
halogen  light  source  and  a  scanning  monochrometer.  A 
particle's  transmission  spectra  and  a  transmission  spectra  of 
a  particle-free  area  of  the  gelatin  coated  slide  are  measured 
and  the  absorption  efficiency  factor,  Qa(A,),  is  calculated  as 

the  difference  normalized  to  the  blank  measurement.  The 
system  is  capable  of  targeting  individual  microscopic 
particles  with  a  lower  size  limit  of  ~3  p,m  spherical  diameter 
(additional  details  are  in  Iturriaga  and  Siegel,  1989). 
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Particulate  Absorption  Coefficients.  di.j^(X)  :  Meltwater  from 
sea  ice  was  filtered  through  glass  fiber  filters  (GF/f ) . 
Filters  were  frozen  and  transported  to  USC  where  spectral 
absorption  coefficients  of  the  bulk  particulates  (ap(?i))  were 

measured  in  a  double  beam  spectrophotometer  (Uvikon-860, 
Kontron  Instr.)  utilizing  a  filter  pad  method  (Truper  and 
Yentsch,  1967;  Mitchell  and  Keifer,  1988).  Scans  of  the 
optical  densities  of  sample  filters  (ODg)  were  run  from  400 
to  7  50  nm  (resolution,  1  nm)  against  a  blank  GF/f  (ODj^)  and 
the  absorption  coefficients  were  calculated  as 


Sp  ( A, )  = 


In (10) (ODs  -  ODb)  A 


A  is  the  clearance  are  of  the  filter,  Vf  is  the  volume  of 
sample  filtered,  and  p  is  the  pathlength  amplification  factor 
(Butler  1962)  resulting  from  scattering  within  the  filters,  p 
was  modeled  according  to  Bricaud  and  Stramski  (1990). 


3n  irradianc 


msmission. 


To  evaluate  the  particulate's  effects 


on  the  apparent  optical  properties  of  the  ice  it  is  necessary 
to  extrapolate  the  measured  inherent  absorption  properties  of 
the  particulates  (Qa's  and  ap's)  to  apparent  absorption 

properties.  Vertically  averaged  absorption  coefficients 
<ap{'k)>  can  be  calculated  as. 


(5) 
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<ap iX) > 


n  * 
Z 


where  ap(^)n  (m“l)  is  the  particulate  absorption  coefficient 
and  Zn  (m)  is  the  thickness  of  each  section  (n)  within  an  ice 
core.  Dividing  <ap(X,)>  by  the  mean  cosine  of  the  downwelling 
irradiance  stream  approximates  a  vertically  averaged 

particulate  attenuation  coefficient  <Kdp(X)>.  is  assumed 

to  be  equal  to  0.65  to  approximate  a  diffuse  radiance 
distribution  within  sea  ice (Grenfell ,  1983;  Arrigo  et  al, 
1992).  <Kdp(^)>  estimated  in  this  manner  represents  a 

weighted  average  over  the  entire  ice  thickness  and, 
therefore,  does  not  rely  on  or  account  for  vertical 
variations  in  particulate  concentrations.  This  approach  also 
neglects  the  scattering  effects  of  the  particles  and, 
therefore,  is  likely  to  be  a  conservative  estimate  of  the 
particles  effects  on  the  ice's  apparent  optical  properties. 

The  influence  of  the  sea  ice  particulates  on  the 
transmission  of  irradiances  into  water  column  can  be 
calculated  and  illustrated  by  subtracting  <Kdp(^)>  from  the 

measured  sea  ice  attenuation  coefficients  and  predicting 
transmitted  irradiances  for  a  "particle  free"  scenario. 
Converting  photon  fluxes  (units,  HE  m"2  to  energy  fluxes 

(units,  W  m“2)  according  to 


W  m-2  =  he  m-2  s-1  * 


(6) 
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{h  =  Planck's  constant,  c  =  speed  of  light)  and  integrating 
the  difference  between  the  measured  spectra  and  the  predicted 
particle-free  scenario,  produces  an  estimate  of  the  reduction 
of  transmitted  energy  due  to  the  presence  of  particulate 
impurities.  A  similar  analysis  can  be  made  for  first-year 
ice  in  the  refrozen  lead.  However,  because  we  were  not  able 
to  measure  attenuation  coefficients  of  the  lead  ice,  we 
utilize  Perovich's  (1990)  approximations  for  first-year  blue 
ice  and  predict  transmission  spectra  with  and  without  the 
addition  of  <Kdp(X)>. 


Results ; 

General  observations:  Salinities  (measured  with  a  Bausch  and 
Lomb  refractometer)  of  the  MY  ice  ranged  from  0  ppt  near  the 
surface  to  15  ppt  near  the  ice  water  interface  -profiles 
indicative  of  MY  ice  undergoing  brine  drainage.  FY  ice 
salinities  ranged  from  5  to  10  ppt  and  exhibited  a  typical 
"C"  shaped  profile  for  FY  ice.  MY  ice  cores  contained 
internal  bands  and  isolated  pockets  of  visibly  detectable 
particulates  while  FY  ice  appeared  particulate-free  to 
unaided  eyes. 

Individual  particulates  and  their  absorption  properties: 
Microscopic  observations  of  particulates  collected  from  both 
the  MY  and  FY  ice  revealed  a  diverse  assemblage  of 
particulates  including:  empty  diatom  frustules,  various 
species  of  pigmented  pennate  diatoms,  tetrads  of  cells  of 


unknown  taxonomic  affinity,  arthropod  remains,  flocculant  9 

detrital  material,  minerals,  and  aggregates  containing  all  of 
the  above  mentioned  material. 

The  majority  (61%)  of  the  algal  cells  were  pennate  diatoms 
ranging  in  size  from  10  to  80  ^im  in  length.  Their  Qa(A.) 

spectra  were  typical  for  marine  diatoms-  having  primary  and 
secondary  absorption  peaks  centered  at  435  and  670  nm  (figure 
lA) .  Within  the  cold  and  low-saline  upper  regions  of  the  MY 
ice,  spherical  spores  (8  to  20  |J.m  diameter),  often  found  in 
tetrads,  were  abundant  among  an  assemblage  of  empty  pennate 
diatom  frustules.  Spores  had  distinct  absorption  spectra 
with  a  very  broad  primary  absorption  peak  centered  at  476  to 
480  |xm  (figure  lA)  -indicative  of  a  large  compliment  of 
photoprotective  and/or  accessory  pigments.  Flagellated 
microalgae  also  were  present  and  were  relatively  more 
abundant  near  the  bottom  of  the  FY  ice  but  were  also  found 
within  the  visibly  discernable  bands  in  the  MY  ice. 

Flagellates  had  primary  absorption  peaks  centered  at  435  and 
665  }im  with  a  pronounced  shoulder  at  470  p.m  (figure  lA) 

-spectral  characteristics  indicative  of  Chlorophyll  b 
absorption  (Bidigare  et  al,  1988) 

The  majority  (70%)  of  the  non-cellular  particulates  were 
grayish  appearing  amorphous  aggregates  of  flocculant 
material.  Qa(^)  were  larger  at  shorter  wavelengths  and  had  a 

monotonic  decrease  towards  the  red  wavelengths  (Figure  IB) . 

This  spectral  signature  is  characteristic  for  marine  detritus 
(Kishino  et  al,  1984;  Roesler  et  al,  1989;  Bricaud  and 
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Stramski,  1990) .  The  remaining  30%  of  the  non-cellular 
particles  were  either  mineral  particulates  (identified  by 
angular  edges,  facets,  color,  and  cross  polarization 
microscopy)  or  very  dense  orange  appearing  amorphous 
particles.  The  orange  appearing  amorphous  particulates  had 
very  high  and  relatively  constant  absorption  efficiencies  at 
wavelengths  less  than  500-515  nm  and  at  wavelengths  greater 
than  500  nm,  a  monotonic  exponential  decrease  in  Qa(X) 

(figure  IB) .  They  appeared  as  dense  orange  particles  when 
viewed  with  white  light  microscopy  but  they  exhibited  no 
other  distinguishing  features.  These  particles  will  be 
designated  as  "type  ii  detritus"  in  this  paper  in  order  to 
distinguish  them  from  the  grayish  detrital  particles  that 
will  be  designated  as  "type  i  detritus".  Several  of  the 
mineral  particulates  (5  to  30  |i.m)  had  distinct  Qq^{X)  spectra 

having  a  concave  upward  signatures  (figure  IB).  Using  cross 
polarizing  microscopy,  several  of  these  minerals  were 
identified  as  biotite  particles.  Feldspar,  quartz,  and  other 
transparent  clay  minerals  were  also  present,  but  their 
absorption  efficiencies  were  below  levels  of  detection. 

Small,  3  -  20  ^im  particles  of  biotite,  chlorite, 
quartz/feldspar,  and  amorphous  soot  particulates  were  also 
abundant  in  snow  samples  examined  via  cross  polarization 
microscopy.  Absorption  efficiency  spectra  were  larger  at 
shorter  wavelengths  with  no  other  apparent  distinguishing 
features  (Figure  1C) . 


Bulk  absorption  coefficients;  Vertical  profiles  of  the 
absorption  coefficients  of  the  bulk  particulates  collected 
from  the  ice  cores  show  that  the  amorphous  type  i  detritus 
was  optically  predominant  throughout  much  of  the  MY  and  FY 
ice  (figure  2  A  and  B) .  Sections  of  ice  that  contained 
visibly  apparent  discolorations  had  increases  in  the  relative 
contribution  of  algal  pigments  to  the  bulk  absorption 
properties  of  the  particulate  matter  (detected  as  a 
discernable  peak  at  675  nm) .  Ap(X)s  of  bulk  samples 

extracted  from  visible  “pockets"  of  particulates  were  much 
larger  than  those  within  the  ice  cores  and  were  dominated  by 
algal  pigments  (figure  2C) .  A  detrital  component  was  still 
apparent  in  these  samples  -  causing  the  primary  absorption 
peak  to  be  centered  at  410  to  420  nm  rather  than  the  435  nm 
chlorophyll  a  peak  measured  for  individual  algal  cells. 

Albedo ;  Albedos  above  CEAREX's  snow  cover  (figure  3)  were 
measured  on  April  9th  under  clear  sky  conditions.  One  source 
of  variability  is  attributable  to  the  effects  of  differing 
solar  angels  during  the  12 -hour  measurement  period.  The 
notable  features  in  the  albedo  are:  1)  they  are  low  (0.4-0. 5) 
over  all  wavelengths  and  2)  albedos  are  larger  at  the  longer 
wavelengths.  These  features  of  the  albedo  are  consistent 
with  particulate  contaminations. 


Tee  and  Snow  Attenuation  Coefficients;  A  K(jsnw(^)  spectrum  12 


is  illustrated  in  figure  4  along  with  the  under-ice 
irradiance  spectra  used  in  its  calculation.  The  removal  of 
the  snow  cover  increased  the  transmitted  energy  by  four- fold. 
The  Kdsnw(^)  values  are  low  compared  to  dry  compact  snow 

(Perovich,  1990)  but  may  be  reasonable  for  an  ice-snow 
mixture  which  was  observed  at  this  site.  Lower  values  may 
also  be  partly  artifactual;  because  specular  reflection  from 
the  air/ice  interface  is  not  accounted  for  in  equation  #1. 

Attenuation  coefficients  for  the  ice  at  both  of  the  MY 
sites  were  similar  in  magnitude  and  spectral  shape  (figure 
5),  and  fell  between  those  previously  reported  for  "white" 
and  "blue"  ice  in  the  Arctic  (Grenfell  &  Maykut,  1977). 

These  coefficients  are  consistent  with  the  visual 
observations  of  the  ice  cores  being  comprised  of  a  "white 
ice"  layer  ca.  10  cm  deep  underlain  by  more  visually 
transparent  ice. 

PRTtAr.lp.' P.  effects  on  transmitted  irradiances;  Particulate 
matter  at  CEAREX  "0"  camp  was  estimated  to  have  had  a 
relatively  small  influence  on  photon  and  energy  flux  into  the 
water  column  (Figure  7  A-C) .  The  estimated  effect  was  a 
reduction  of  transmitted  energy  by  20%  (0.2  W  m"2  when 
incident  flux  is  300  W  m”2)  through  MY  ice  and  10%  (3  W  m-2) 
through  the  FY  ice.  The  highly  concentrated  pockets  of 
material  were  estimated  to  cause  the  largest  alterations  of 
the  transmitted  energy  and  irradiance  spectra  (figure  7C) . 


The  estimated  effect  was  to  reduce  the  transmitted  energy  by  13 
65%  (5.5  W  m“2)  and  shift  the  wavelength  maxima  from  496  nm  to 
520-526  nm.  This  is  a  comparable  effect  as  those  measured 
and  modeled  for  a  microalgal  blooms  in  Arctic  bottom  ice  near 
the  seasonal  maxima  (Perovich  et  al,  1993). 

Discussion ; 

Albedo ; 

The  snow's  albedo  at  "O"  camp  was  characteristic  of  a  snow 
pack  contaminated  with  particulate  matter  (figure  3) .  Our 
microscopic  analysis  showed  the  presence  of  several  types  of 
amorphous  flocculant  and  lithogenic  particles.  The 
possibility  that  some  of  these  particulates  originated  from 
within  the  camp  cannot  be  entirely  discounted.  However, 
several  investigators  have  found  black  carbon,  minerals,  and 
other  particulates  within  Arctic  snow  and  within  the  famous 
"Arctic  Haze".  Our  spectral  measurements  of  their  absorbing 
properties  (fig  4c)  are  consistent  with  those  measured  at 
discrete  wavelengths  (435,  525,  660  and  880  nm)  on  different 
particulate  assemblages  found  throughout  the  Arctic  Basin 
(Clark  and  Noone,  1985). 

irradiance  Transmission; 

As  has  been  demonstrated  on  multiple  occasions,  snow  cover 
is  a  major  constituent  influencing  the  transmission  of 
irradiances  into  and  through  the  ice  cover  (Grenfell  and 
Maykut,  1975;  Perovich,  1990).  The  snow  cover  at  CEAREX  "0" 
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camp  was  no  exception.  The  snow  cover  accounted  for  an  70- 
90%  reduction  of  shortwave  radiation  (400-700nm)  passing 
through  MY  ice  and  accounted  for  an  estimated  50-80% 
reduction  of  energy  passing  through  FY  ice  in  a  refrozen  lead 
(data  not  shown) . 

Particulate  impurities  in  the  ice  also  contributed  to  the 
reduction  of  transmitted  irradiances  (figure  7  A-C) .  Their 
effect  was  much  less  dramatic  than  that  of  the  snow  cover. 
However,  the  magnitude  of  the  particle's  effect  may  be 
substantial  in  terms  of  energy  transfer  processes  under 
certain  conditions.  For  example,  the  particle's  reduction  of 
short-wave  energy  flux  through  the  lead  under  clear  skies  is 
of  the  same  magnitude  as  the  mean  oceanic  heat  flux  for  the 
central  Arctic  Ocean  (2  W  m"^,  Maykut  and  Perovich,  1987). 

Absorption  efficiency  spectra  of  individual  particulates 
can  be  utilized  in  radiative  transfer  models  to  elucidate  the 
potential  effects  different  types  of  particles  can  on  energy 
transfer  processes.  The  absorption  coefficients  of  the 
particulate  matter  can  be  calculated  from  the  Qa(^) 's  as 

follows : 


=  Q^{X)  *  N  *  g 

where  N  is  the  concentration  of  particles,  and  g  is  their 
optical  cross  section  (Iturriaga  and  Siegel,  1989). 

Different  compliments  of  particulate  loads  are  predicted  to 
change  the  spectral  quality  of  transmitted  irradiances.  For 
example,  the  difference  between  sea  ice  containing  diatoms 
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versus  flagellates  at  equal  cell  numbers  is  a  definitive 
spectral  shifting  of  the  wavelength  maximum  of  the 
transmitted  downwelling  irradiance  spectrum  (fig  8a).  A 
detrital-dominated  particulate  compliment,  on  the  other  hand, 
produces  a  transmitted  downwelling  irradiance  spectrum  with  a 
broad  peak  but  little  spectral  shifting  (fig  8b). 

As  algorithms  are  being  developed  to  remotely  sense  sea  ice 
and  snow  thicknesses  as  well  as  algal  standing  crops  (e.g. 
Legendre  and  Gosselin,  1991;  Perovich  et  al,  1993), 

spectra  may  be  useful  in  the  algorithm  development  for  the 
deconvolution  of  the  bulk  apparent  optical  properties  into 
the  respective  inherent  optical  properties  of  the  sea-ice 
constituents . 

Particles  effects  on  Primary  Production; 

Sea  ice  particulates,  in  effect,  compete  with 
photosynthetic  organisms  for  photons.  Photons  absorbed  by 
ice,  water,  soot,  biotite,  detritus  and  other  non¬ 
photosynthetic  constituents  are  no  longer  available  for 
absorption  by  algae.  Therefore,  the  initial  effects  of 
impurities  in  ice  are  to  reduce  the  potential  rates  of 
primary  production.  In  sea  ice  ecosystems  and  in  the 
underlying  water  column  light  absorption  is  considered  to  be 
a  rate-limiting  ecological  process  (Gross!  et  al.,  1987). 
However,  if  sea  ice  particulates  accelerate  ice  decay  by 
contributing  to  a  positive  feed-back  process  involving  ice 
deteriation  (Eicken  et  al.,  1991)  and  decreasing  ice 
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concentrations  leading  acclerated  ice  decay  (Maykut  and 
Perovich,  1987),  then  water  column  blooms  (and  primary 
production)  may  be  enhanced.  The  potential  positive  and 
negative  feebacks  process  associated  with  these  scenarios 
make  it  difficult  to  predict  the  effects  that  particulates 
have  on  annual  production  cycles.  However,  they  undoubtedly 
play  a  role  in  the  timing  and  location  of  production. 

Energy  Rudaets  and  Regional  Albedos;  The  absorption  of 
electromagnetic  energy  by  particulates  on  the  surface  of  the 
ice  and  snow  has  been  recognized  as  an  important  influence  on 
albedo  (this  study,  Wiscombe  and  Warren,  1982)  and  localized 
energy  budgets.  If  particulates  material,  both  lithogenic 
and  biological  are  present  in  the  ice  over  large  scales  (and 
there  is  evidence  that  thery  are)  then  they  are  potentially 
important  for  local,  regional,  and,  perhaps,  global  energy 
budgets . 

The  potential  for  discrete  bands  and  pockets  of  absorbing 
particulates  within  the  interior  sea-ice  matrix  for  altering 
localized  energy  transfer  processes  and  regional  energy 
budgets  through  indirect  means  is  recognized.  Internal  bands 
of  particulate  matter  leading  to  internal  heating  or  melting 
will,  in  turn,  influence  the  optical  (Grenfell,  1983)  and 
mechanical  (Buynitsky,  1968)  properties  of  ice.  Therefore, 
it  is  likely  that  internal  heating  contributes  to  the 
positive  feedback  dynamic  that  accelerates  the  decay  of  a 


sea-ice  cover. 


Accelerated  opening  of  ice-free  water,  in 
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turn,  causes  lowered  regional  albedos. 

Algal  communities  on  the  bottom  of  FY  Arctic  ice  as  well  as 
Antarctic  FY  ice  have  been  observed  to  ablate  earlier  than 
the  surrounding  ice  due  to  localized  heating  effects  of  the 
algae  (Gross!  et  al,  1987;  Lewis  and  Cota,  unpublished). 
Cryoconite  holes  forming  on  the  surface  of  the  ice  are 
another  example  of  particulates  having  an  easily 
recogni zable / ident i f iable  influence  on  ice  dynamics. 

Internal  heating  and  melting  within  sea  ice,  however,  will 
only  manifest  its  effects  as  alterations  of  internal  sea  ice 
structure  and/or  accelerations  of  ice  decay  at  the  lower 
boundary  and,  therefore,  is  a  phenomena  that  is  less  readily 
observable . 

Multiple  studies  indicate  particle-laden  pack  ice  in  the 
Arctic  basin  is  common  and  is  a  significant  carrier  for 
cross-shelf  transport  of  sediments  and  continental  shelf 
erosion  (Pfirman  et  al,  1989).  The  potential  significance  of 
ic  le  -  laden  ice  throughout  the  Arctic  basin  on  energy 
budgets  has  been  suggested  (e.g.  Pfirman  et  al.,  1987)  and 
SiXQ.  particularly  reinforced  by  our  findings  of  internal 
banding  and  their  potential  direct  and  indirect  effects  of 
internal  heating.  Because  particulate  material  in  sea  ice 
and  snow  is  variable  in  both  space  and  time  and  the  pack  ice 
is  relatively  inaccessible,  we  can  not  yet  determine- 
guant itatively—  the  role  internalized  particulates  have  in 
determining  the  optical,  thermal,  and  biological  properties 


of  Arctic  sea  ice  on  regional  scales.  Further  analysis  of  the  18 
factors  controlling  the  spatial  (both  vertical  and 
horizontal)  and  temporal  distribution  of  terrigenous  matter, 
aeolian  particulates,  detritus,  and  algae  in  sea  ice  seems 
warranted.  Moreover,  the  development  of  remote  sensing 
algorithms  for  detecting  particulate  loads  in  ice  will  be 
useful  tools  in  understanding  the  dynamics  of  polar 
ecosystems.  The  presence  of  several  different  types  of 
particulates  with  differing  optical  characteristics  makes  the 
interpretation  of  existing  algorithms  (Legendre  and  Gosselin, 
1991;  Perovich  et  al,  1993)  problematic. 


Figure  Legends 
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Figure  1.  Absorption  efficiency  spectra  of:  A),  diatoms, 
flagellates  and  spores  found  within  MY  pack  ice  (Flagellates 
and  diatoms  were  also  common  in  FY  ice  samples) ,  B) . 
detrital  and  mineral  particulates  found  within  MY  pack  ice, 
and  C) .  snow  particulates. 

Figure  2.  Bulk  particulate’s  absorption  coefficients  of 
particles  collected  from  discrete  sections  of:  A)  MY  ice 
cores;  B)  FY  ice  cores;  and  C)  visibly  apparent 
discolorations  within  MY  ice  -illustrated  as  bands  and 
pockets  of  dark  matter. 

Figure  3.  Spectral  Albedos  of  Snow.  Measured  values 
(mean  ±  2  standard  errors,  n=  8)  are  those  for  a  X  cm 
metamorphosed  snow  cover  (i.e.  firn)  at  CEAREX  "0"  camp. 

Pure  snow,  100  ppmw  ash,  and  1,000  ppmw  ash  spectra  are 
albedos  predicted  by  Wiscombe  and  Warren  (1982)  for  varying 
concentrations  (ppm  by  weight,  ppmw)  of  Mt.  St.  Helen's  ash 
in  a  model  snow  pack  with  grain  radii  of  1,000  ^.m. 

Figure  4.  A).  Downwelling  spectral  irradiances  beneath 

2.x  meters  of  MY  ice  with  and  without  20  cm's  of  snow. 
Incident  surface  irradiances  were  1710-1720  p.E  m“2  and 

the  downwelling  irradiance  beneath  the  ice  was  0.3  p.E  m~2 
(0.073  W  m“2)  with  snow  and  1.39  JXE  m“2  s-1  (0.329  W  m“2) 


without.  B)  .  Diffuse  attenuation  coefficient  for  snow  (n 
13,  Error  bars,  +  2  standard  deviations). 

Figure  5.  Diffuse  attenuation  coefficients  for  MY  ice  at 
"0"  camp.  Also  shown  are  attenuation  coefficients  for  "blue" 
ice  and  "white"  ice  (after  Grenfell  and  Maykut,  1975)  . 

Figure  6.  Predicted  downwelling  spectral  irradiances 
beneath:  A) .  MY  ice  with  a  20  cm  snow  cover  with  and  without 
<Kp>s  calculated  for  "0"  camp  MY  ice.  Predictions  suggest 
that  particulates  accounted  for  a  20%  reduction  (0.5  W  m"2 
less  when  incident  irradiance  =  300  W  m"2)  of  transmitted 
energy;  B) .  FY  ice  with  X  cm  snow  cover  with  and  without 
<Kp>s  of  "0"  camp  FY  ice.  Particulates  account  for  a 
reduction  of  3.0  W  m"2  -  a  10%  reduction  of  transmitted 
energy;  C) .  2.78  m  of  MY  ice  with  a  10  cm  band  of  particles 

with  aps  measured  in  algal  "pockets"  (fig.  6c) .  Energy 
reduction  due  to  a  band  of  particulates  and  algae  =  5.5  W 
-  a  65%  reduction  in  energy  transmission. 

Figure  7 .  Predicted  downwelling  spectral  irradiances 
beneath  ice  containing  comparable  numbers  of  A) .  diatoms  or 
flagellates;  B) .  diatoms  or  detritus. 
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ABSTRACT 

Time  series  studies  of  the  spectral  irradiance  fields  beneath  multiyear  pack  ice  were  conducted  in  the  Eastern  Arctic 
Basin  at  82-8  30N  as  part  of  the  Coordinated  Eastern  Arctic  Research  Experiment  (CEAREX).  Particulate  matter  was 
collected  from  the  multiyear  pack  ice  as  well  as  first  year  ice  in  a  refrozen  lead.  The  vertical  distribution  within  the  ice  and 
the  spectral  absorbtion  properties  of  the  particulates  were  determined  in  order  to  estimate  their  contribution  to  the  optical 
properties  of  the  sea  ice. 

Among  the  particulates  contained  in  sea  ice  detritus  was  common  throughout  all  portions  of  the  pack  ice  and  was 
the  major  light  absorbing  particulate  matter  in  the  ice  at  the  time  of  the  observation.  Algal  cells  and  mineral-like 
particulates  also  were  present,  yet  they  contributed  to  the  light-absorbing  properties  to  a  lesser  extent  than  the  detritus. 
During  early  spring,  particulate  matter  contributed  litde  to  the  bulk  attenuation  coefficients  of  the  multiyear  ice,  however,  it 
was  estimated  to  have  a  more  substantial  contribution  to  the  attenuation  coefficients  of  first  year  ice  in  a  refrozen  lead. 
Results  of  a  single  stream  multilayer  radiative  transfer  model  that  simulates  concentrations  of  biogenic  particulate  matter 
observed  in  Arctic  sea  ice  indicates  that  particulate  matter  within  sea  ice  plays  a  substantial  role  in  radiative  energy  transfer 
and  has  the  potential  to  seasonally  alter  spectral  irradiance  regimes  within  the  ice  covered  Arctic  ocean. 

1.  Introduction 

In  the  ice  covered  regions  of  the  polar  oceans,  the  transfer  of  electromagnetic  radiation  (EMR)  into  the  upper  layers 
of  the  water  column  is  strongly  dependent  on  the  inherent  optical  properties  of  the  sea  ice  and  snow.  Particles,  which  are 
sometimes  referred  to  as  contaminants  L2  will  conuibute  to  the  inherent  optical  properties  of  the  ice  cover  to  varying 
degrees  depending  on  their  concentrations  and  their  absorption  and  scattering  properties.  The  influence  of  particles  on  the 
transmission  of  irradiance  through  sea  ice  has  been  observed  in  the  field  by  many  woikers.3.4.5.6 

Recently,  two  modeling  efforts  have  incorporated  the  effects  which  particulate  matter  has  on  the  propagation  of 
EMR  in  sea  ice.  Arrigo  et  al7  developed  and  then  tested  with  field  data  a  bio-optical  model  describing  seasonal  changes  in 
downwelling  spectral  inadiances  (400-700  nm)  within  and  beneath  land-fast  congelation/platelet  ice  containing  biogenic 
particles,  (jrenfell*  included  absorption  and  scattering  terms  for  impurities  in  a  two-stream  multi-layer  radiative  transfer 
model  for  sea  ice.  These  models  indicate  that  absorption  properties  of  different  species,  taxonomic  groups,  or  assemblages 
of  microalgae  may  cause  substantial  variability  in  the  quantity  and  quahty  of  under-ice  spectral  inadiances  within  the  visible 
spectrum  (4(X)-700  nm). 

During  the  1989  Coordinated  Eastern  Arctic  Experiment  (CEAREX)  we  collected  particulate  matter  from  the  sea 
ice  North  of  Svalbard  Norway  during  the  early  spring.  TTie  relative  abundance  and  morphology  of  particle  types  and  their 
light  absorbing  properties  were  determined.  A  concurrent  assessment  was  made  of  ie  ambient  surface  and  under-ice 
spectral  irradiance  fields  in  order  to  estimate  the  influence  of  the  particulate  matter  on  the  apparent  optical  properties  of  the 


sea  ice.  Single  stream  multilayer  radiative  transfer  models  are  utilized  to  estimate  the  influence  these  particles  play  in 
energy  transfer  and  alterations  of  spectral  irradiance  fields. 


2.  METHOn.S 


2.1.  Irradiance  measurements 

During  operations  at  the  CEAREX  oceanographic  camp  (designated  "0"  camp)  (originally  located  at  820n,  10®W) 
irradiance  measurements  were  made  from  April  7  to  AprU  22.  Measurements  included:  under-ice  spectral  downweUing 
irradiance  (E(j(z,A.)),  surface  downweUing  spectral  irradiance  (Ed(0+A),  and  surface  downweUing  scalar  irradiance 
(Eod(0-t-,PAR))  as  photosynthetically  available  radiation  (PAR)  (4(X)-700nm).  Eod((>+,PAR)  was  measured  with  a 
Biospherical  Instruments  model  QSR-2(X)  Solar  Reference  Hemispherical  Irradiance  sensor.  Spectral  measurements  were 
made  with  a  MERlOlO  spectroradiometer  (Biospherical  Instruments)  with  its  5  nm  spectral  band  widths  centered  at  410, 
441  455,  488,  507,  532,  550,  570,  589,  633,  656,  and  694  nm.  Both  instruments  were  calibrated  before  and  foUowing  the 
CEAREX  89  experiments  at  Biospherical  Insmunent  company  in  San  Diego,  California. 

Under-ice  measurements  with  the  undisturbed  snow  cover  intact  were  made  by  mooring  the  spectroradiometer  on  an 
articulated  arm  deployed  through  a  75  cm  diameter  hydrohole  (HH)  drilled  with  a  steam  heat  drill.  The  instrument  was 
positioned  at  the  ice/water  interface  approximately  3  m  away  from  the  hole  which  was  covered  to  avoid  stray  light  from 
reaching  the  cosine  coUector.  With  the  radiometer  still  in  position,  natural  snow  cover  was  manually  removed  from  above 
the  radiometer  in  order  to  estimate  the  bulk  diffuse  spectral  attenuation  coefficients  for  the  sea  ice  and  snow  using  methods 
described  below. 

2.2.  Sample  collections 

Two  study  sites  were  established  where  under-ice  inadiances  were  measured  and  ice  cores  were  also  obtained  These 
will  be  referred  to  as  HHl  and  HH2.  HHI  and  HH2  were  located  in  a  relatively  flat  area  among  hummocked  ice.  Ice 
thickness  was  2.75  meters  and  snow  cover  was  20  to  30  cm.  Average  salinities  were  3-4  ppt  (data  not  shown).  Ice 
thickness,  surface  topography,  snow  cover,  and  low  salinities  indicate  this  was  multiyear  ice.  An  additional  study  site  was 
established  in  an  undeformed  refrozen  lead  (L)  with  1.2  meter  thick  ice  having  10  to  20  cm  snow  cover.  Two  ice  cores, 
designated  LI  and  L2  were  obtained  from  this  site  with  average  salinities  of  5-6ppt.  The  relatively  high  salinities,  lack  of 
deformation,  and  ice  and  snow  thickness  indicate  that  this  was  ice  which  had  formed  in  the  previous  winter  (i.e.  first  year 
ice). 


Ice  cores  were  obtained  using  a  CRREL  auger  (7  cm  diameter),  sectioned,  and  then  allowed  to  melt  at  2^C  in 
closed  opaque  containers.  Samples  for  the  analysis  of  the  absorption  properties  of  the  particulate  material  were  obtained 
from  this  meltwater  and  were  analyzed  according  to  methods  described  below. 

2.3.  Diffuse  Attenuation  Coefficient  for  Snow.  KdsnwpLi 


When  the  snow  cover  was  removed  from  above  the  radiometer,  the  diffuse  spectral  attenuation  coefficient  for  the 
natural  snow  cover  K<jsnw(^)  (m-l)  was  estimated  using  the  Beer-Lambert  model  describing  the  exponential  attenuation  of 
hght  through  a  homogenous  medium. 


Kdsnw(^)  ~  "h  *  ^ 


E<i(-t-snw,A,) 

E(j(-snw,X) 


(1) 


where  H  is  the  thickness  of  the  snow  in  meters  and  E(j(+snw,X,)  and  E(j(-snw,A.)  are  the  under-ice  downweUing  spectral 
irradiances  ()iE  m'^  s'^)  with  and  without  snow,  respectively. 

2.4,  Diffuse  attenuation  coefficient  for  sea  ice.  K^jpp(X) 

Also  using  the  Beer-Lambert  model,  the  bulk  diffuse  spectral  attenuation  coefficient  for  optically  thick  sea  ice 
(K<iice  (^))  can  be  approximated  from  the  following  relationship; 


(2) 


.  E<j(z.X) 
E^(O-A) 


(9,10)  where  z  is  the  thickness  of  the  ice  in  meters,  E(j{z,X)  is  the  downwelling  irradiance  at  the  bottom  of  the  ice  and  £^(0- 
X)  is  the  downwelling  irradiance  just  below  the  surface  of  the  ice.  EdfO-.X)  cannot  be  measured  direcUy  in  the  field  but  can 
be  estimated  as 


Ed(0-,>.)  =  EdCO+.X)  *(l-Ro) 

where  Rq  is  the  specular  reflection  from  the  sea  ice  surface  which  is  estimated  accca-ding  to  Arrigo  et  al.7 
2.5.  Attenuation  coefficient  of  particulate  matter. 


Attenuation  of  downwelling  irradiance  in  a  medium  is  due  to  the  combined  attenuation  due  to  the  medium  itself 
and  attenuation  by  particles  and  dissolved  substances  within  the  medium. 1E12  Ignoring  attenuation  by  dissolved 
substances,  the  diffuse  spectral  attenuation  coefficient  for  sea  ice  may  be  approximated  by: 

Kdice(^)  =  Kdi(A.)  +  Kdp(A.)  ^4] 

where  Kdp(X)  is  the  spectral  diffuse  downwelling  attenuation  coefficient  of  particulate  material  and  Kdi(X)  is  the  spectral 
diffuse  downwelling  attenutation  coefficient  of  pure  sea  ice.  Furthermore,  Kdp(X)  can  be  approximated  by 


Kdp(^)  = 


^p(^)  +  bp(A.) 
ud(z) 


(5) 


(13)  where  ap(X)is  the  total  particulate  absorption  coefficient  (m-l),  bp(X)  is  the  particulate  backscattering  coefficient  (m-l) 
and  ud(z)  is  the  mean  cosine  of  the  angular  distribution  of  downwelling  irradiance  which  is  assumed  to  be  0.656.^^ 
Backscattering  by  sea  ice  particulate  matter  was  not  measured  and  wiU  not  be  taken  into  account  at  the  present  time  in  order 
to  approximate  K<ip(X) . 

16^  Particulate  Absorption  Coefficient  a^^tX) 

In  order  to  estimate  spectral  light  absorbing  properties  of  the  particulate  matter  contained  within  the  ice,  ice  core 
meltwater  was  filtered  through  glass  fiber  filters  (GF/f)  and  the  spectral  absorption  coefficient  of  the  total  particulate  matter 
(ap(X))  was  estimated  utilizing  filter  pad  methods  which  have  been  described  previously.5.15,16 

2.7,  Absorption  Efficiency  of  Individual  Particles.  O^fX.) 

In  order  to  determine  the  light  absorbing  fH'operties  of  individual  particles,  ice  core  meltwater  was  filtered  through 
0.4  pm  pore  size  Nuclepore  filters  at  low  pressures  (5  mmHg).  Particles  retained  by  the  filters  were  immediately 
transferred  to  gelatin  coated  microscope  slidesl^  and  kept  frozen  until  they  were  analyzed  at  USC.  Absorption  efficiency 
spectra  of  the  individual  particles  were  determined  using  a  microphotometric  technique.  1^ 

2.8  Particule's  contribution  to 


A  vertically  averaged  absorption  coefficient  <ap(A.)>  of  the  sea  ice  particulate  matter  can  be  calculated  as, 

2^p(^)n*  Zn 
Z 


(6) 


where  ap(X)n  (m-D  is  the  particulate  absorption  and  Zn  (m)  is  the  thickness  of  each  section  (n)  of  the  ice  core.  By 
substituting  <ap(X.)>  for  ap(X)  in  equation  10  a  vertically  averaged  particulate  attenuation  coefficient  <K<ip(A,)>  for  the  ice 
cover  can  be  approximated.  <K<jp(X)>  estimated  in  this  manner  represents  a  weighted  average  over  the  entire  ice  thickness 


and.  therefore,  does  not  rely  on  or  account  for  vertical  variations  in  particulate  concentrations.  <K<jp(X)>  can  then  be  used 
to  approximate  the  particulate  matter's  contribution  to  the  sea  ice  attenuation  coefficients,  F(?t),by 


Fa)  = 


<K<jf)(A.)> 

K<jice(^) 


2.9.  Modeling  energy  absorption  within  sea  ice 


(7) 


EMR  absorption  by  the  constituents  of  a  homogeneous  medium  depends  on  the  constituent's  absorbing  properties, 
their  vertical  distribution,  and  the  irradiance  distribution  within  the  medium.  Given  the  absorption  coefficient  of  constituent 
X,  the  attenuation  coefficient  (K<i(A.))  for  the  medium  and  assuming  a  homogenous  distribution  of  each  constituent  within 

the  medium,  the  rale  of  EMR  absorption  by  constituent  x  fE^^(A.))  can  be  approximated  by 


E^\s(^)=[Ed(Oa)-Ed(Z.X)]*^  (8) 

where  E(i(0-,X)  is  the  downwelling  irradiance  just  below  the  surface  of  the  medium,  Ed(Z,A.)  is  the  transmitted  downwelling 
spectral  irradiance.  The  ratio  of  the  absorption  coefficient  to  the  attenuation  coefficient  simply  is  the  fraction  of  attenuated 
downwelling  irradiances  lost  to  absorption  by  constituent  x.  Ed(Z,X)  can  be  determined  according  to  the  Beer-Lambert 
model. 


Ed(0-,X)  *  e(-K<i(^)  *  Z).  (9) 

using  a  finite  integration  scheme 
700 

=  /  ^abs^^^*^  (10) 

h=ix) 

in  order  to  obtain  the  rate  of  absorption  of  total  visible  EMR  by  each  constituent  within  the  medium.  Because  sea  ice 
usually  consists  of  multiple  layers  with  distinct  optical  properties  equations  8,9,  and  10  are  applied  for  multiple  layers 
when  modeling  absorption  within  a  sheet  of  sea  ice.  At  the  surface  of  the  ice  or  snow^  specular  reflection  is  modeled 
according  to  Arrigo  et  al^  and  specular  reflection  at  the  interface  between  layers  is  assumed  to  be  zero.  1-7  For  our  purposes, 

X  is  considered  to  be  either  pure  sea  ice  (i)  or  particulate  material  (p).  We  wished  to  model  the  partiUoning  of  visible  EMR 
that  was  absorbed  by  the  sea  ice  and  particul^s  within  the  ice,  and  the  EMR  passing  into  the  water  column.  This  was 
accomplished  by  combining  equations  3,4,  8,  and  9  and  specifying  an  incident  spectral  irradiance  (400-700  nm)  (figure  2), 
attenuation  coefficients  for  snow  (figure  lb),  absorpuon  coefficients  for  the  sea  ice  as  "bubble  free"  sea  ice,  19  attenuation 
coefficients  of  sea  ice,  1-3.7,  and  absorption  coefficients  for  the  particulates  within  the  ice. 


Ed(Z,X)  = 


Eabs  (^)  can  be  spectrally  integrated 


"abs 


3.  RESULTS 

31.  Diffuse  Attenuation  Coefficient  for  Snow,  (X) 

Under-ice  spectral  downwelling  irradiances  on  Julian  day  97  and  98  (April  7  and  8)  with  and  without  the  natural 
snow  cover  intact  are  shown  in  figure  la  Using  these  measurements,  K<isnw  (^)  was  estimated  according  to  equation  #1 
(figure  lb).  The  shape  of  the  resultant  K(isnw(^)  spectrum  is  typical  for  snow^  however,  the  magnitude  is  low  compared 

to  dry  compact  snow^O  but  may  be  reasonable  for  an  ice-snow  mixture  which  was  observed  at  this  site.  The  low  values 
may  also  be  an  artifact  because  specular  reflection  from  the  air/ice  interface  is  not  taken  into  account  in  equation  #1 . 
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Figure  1.  A)  Downwelling  spectral  irradiances  beneath  2.4  meters  of  multiyear  ice  with  and  without  20  centimeters  of 
snow.  Eod(0+,PAR)  =  1721  and  1711  nE  m'2  s'l  during  the  measurement  periods  with  snow  and  without  snow, 
respectively.  B)  Estimated  diffuse  attenuation  coefficient  of  the  snow  cover  based  on  spectra  shown  in  A. 

3.2.  Diftiise  attenuation  coefTicients  for  sea  ice.  Kdice 

Diel  time  series  of  surface  downwelling  spectral  irradiances  on  JD  99  under  clear  skies  are  shown  in  figure  2a. 
When  normalized  to  their  peak  wavelengths  (figure  2b)  it  is  apparent  that  spectral  shifting  towards  the  shorter  wavelengths 
occured  at  increased  solar  angles.  This  is  consistent  with  a  shifting  toward  a  more  diffuse  skyhght  which  is  predominantly 
blue.2i  Under-ice  downwelling  spectral  irradiances  at  the  HHl  position  are  shown  in  figure  3  with  the  snow  cover  removed 
and  the  sky  state  was  similar  to  that  of  JD  98  when  measurements  of  the  surface  downwelling  spectral  irradiances  were 
made.  Note  that  spectral  shifting  towards  the  blue  wavelengths  also  occurs  beneath  the  ice  (figure  3b). 

Kdice(^)  spectra  for  the  HHl  and  HH2  study  sites  are  illustrated  in  figure  4.  These  values  fall  between  those 
reported  for  a  "white  ice"  granular  or  scattering  layer  or  those  for  the  interior  portions  of  multi-year  "white  ice"  or  blue 
ice.22  This  is  consistent  with  visual  observations  made  in  the  field  where  it  was  noted  that  the  top  portions  of  75  cm 
diameter  ice  cores  were  indeed  white  while  the  remainder  of  the  cores  appeared  blue-green.  Therefore,  we  recognize  that  the 
ice  was  not  a  structurally  homogeneous  medium  but  instead  had  at  least  two  distinct  layers  that  could  be  visually 
distinguished. 
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Figure  2.  A)  Diurnal  time  series  of  incident  downwelling  spectral  irradiances  measured 
between  0800  to  2300  hours  on  JD  99  B)  Spectra  were  normalized  to  the  peak  wavelength 
in  order  to  illustrate  spectral  shifting. 
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Figure  3.  A)  Diurnal  time  series  of  downweiling  spectral  irradiances  beneath  2.75  meters 
of  multiyear  pack  ice  (HHl  sampling  site)  recored  on  Julian  day  98.  B)  Spectra  were 
normalized  to  the  peak  wavelength  in  order  to  illustrate  spectral  shifting. 
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Figure  4.  Spectral  diffuse  attenuation  coefficients  for  the  multiyear  ice  at  HHl  (open  squares)  and  the  HH2  (solid 
squares)  sites. 


3.3.  Individual  particle's  absorption  efficiencies.  0^(X) 

The  ice  cores  from  HHl  and  HH2  contained  faint  bands  of  particulate  matter,  one  meter  below  the  ice/snow 
interface,  visible  to  unaided  eyes.  Particulate  matter  was  not  visually  apparent  in  the  remaining  sections  of  these  cores  nor 
in  the  cores  from  the  refrozen  lead.  Microphotometric  analysis  revealed  that  all  sections  of  the  multiyear  and  lead  ice  cores 
contained  a  variety  of  fine  grained  (<500  )im)  particulates.  The  multiyear  ice  contained  a  diverse  assemblage  of  particulates 
including:  empty  diatom  frustules,  various  species  of  pigmented  pennate  diatoms,  tetrads  of  cells  of  unknown  taxonomic 
affinity,  arthropod  remains,  flocculant  detrital  material,  individual  mineral-like  particulates,  and  aggregates  containing  all  of 
the  above  mentioned  material.  The  ice  from  the  refrozen  lead  also  contained  particulates  similar  to  that  found  in  the  HHl 
and  HH2  ice  cores.  Some  of  the  microalgae  found  in  the  bottom  sections,  however,  were  morphologically  distinct  from 
those  found  in  the  HH 1  and  HH2  ice  cores.  Total  pigment  (chlorophyll  a  and  pheopigment)  concentrations  were  less  than 
0.15  mg  m"3  in  all  sections  of  both  the  multiyear  ice  and  ice  from  the  refrozen  lead. 

Qa(>.)  spectra  of  the  tnicroalgae  found  in  the  HHl  ice  core  are  shown  in  figure  5.  Tlie  majority  of  the  cells  (61%) 
were  pennate  diatoms  ranging  in  size  from  10  to  80  lun  in  length.  Their  Qa(X)  revealed  shapes  typical  of  spectra  for 
marine  diatoms  having  primary  and  secondary  absorption  peaks  at  430  and  670  nm  respectively. 23  There  were,  however, 
some  pigmented  cells  with  atypical  absorption  spectra.  Most  noteworthy  were  some  pigmented  spherical  cells  found  in 
section  A  which  were  often  found  in  tetrads.  The  cells  were  8  to  20  |4m  in  diameter  and  their  absorption  spectra  were 
characterized  by  broad  peaks  centered  around  470  to  480  nm.  Furthermore,  a  few  of  the  diatoms  from  the  lower  sections  of 
the  ice  possessed  very  high  absorption  efficiencies  with  broad  and  flat  primary  and  secondary  absorption  peaks  indicative  of 
packag^  pigments  or  large  cells.23 

The  non-algal  particulates  possesed  Q^iX)  spectra  with  a  wide  variety  of  shapes  (figure  5).  The  majority  (70%) 
were  amorphous  aggregates  of  flocculant  material  with  Qa(A.)  being  larger  at  shorter  wavelengths  and  having  a  monotonic 
deaease  towards  the  red  wavelengths  which  is  typical  for  detritus.24  The  remaining  30%  of  the  non-algal  particles  appeared 
to  be  mineral-like  particulates  because  of  their  characteristic  angular  edges,  facets,  and  color.  These  also  had  high 
absorption  efficiencies  at  shorter  wavelengths.  However,  unlike  the  detritus,  their  Qa(A.)  spectra  were  concave  upward  and 
were  relatively  high  and  flat  from  400  to  525  nm  with  a  rapid  dechne  at  the  longer  wavelengths.  Mineral-like  particulates 
ranged  from  5  to  15  (im  in  diameter.  Similar  Qa(A.)  spectra  were  measured  from  individual  particulate  matter  from  HHl  and 
HH2  ice  COTes. 


3.4.  Particulate  matter  absorption  coefficients  ap(X) 

Absorption  coefficients  of  the  total  particulate  matter  collected  from  the  ice  cores  revealed  that  detritus  was  the 
predominant  light  absorbing  particulates  throughout  vertical  profiles  of  both  ice  types  (figures  6,7,8,9).  Section  B  of  the 
HHl  ice  core,  was  the  only  portion  of  this  core  which  had  ap(X)'s  with  a  distinct  chlorophyll  absorption  peak  at  675  nm. 
All  of  the  HH2  sections,  on  the  other  hand,  had  ap(A,)'s  with  distinct  absorption  peaks  at  675  nm.  The  ap(X)'s  from  the  LI 
and  L2  ice  cores  showed  the  675  nm  absorption  peak  only  in  the  bottom  most  sections  of  the  cores.  The  HH2  ice  core  was 
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Figure  5.  Absorption  efficiency  spectra  of  individual  pigmented  cells  (left  column)  and 
non  pigmented  particles  (right  column)  collected  from  discrete  sections  (A  through  D)  of 
the  multiyear  ice  core  (HHl).  Section  locations  within  the  HHl  ice  core  and  thickness  are 
illustrated  in  figure  6. 
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Figure  6.  Diagram  illustrating  the 
multiyear  ice  core  (HHl)  and  the 
particulate  absorption  coefficients  as  a 
function  of  wavelength  (400-700nm)  for 
the  particulates  collected  from  vertical 
sections  A  through  D.  A  5  cm  thick  algal 
layer  is  illustrated  at  1.75  meters  from 
the  ice/snow  intervace  which  was 
visually  apparent  in  the  field. 


Figure  7  Diagram  illustrating  the 
multiyear  ice  core  (HH2)  and  the 
particulate  absorption  coefficients  as  a 
function  of  wavelength  (400-700nm)  for 
the  particulates  collected  from  vertical 
sections  A  through  D.  A  5  cm  thick  algal 
layer  is  illustrated  at  1.75  meters  from 
the  ice/snow  intervace  which  was 
visually  apparent  in  the  field. 
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Figure  8  Diagram  illustrating  an  ice 
core  (LI)  from  a  refrozen  lead  and  the 
particulate  absorption  coefficients  as  a 
function  of  wavelength  (400-700nm)  for 
the  particulates  collected  from  vertical 
sections  A  through  D  of  the  core. 
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Figure  9  Diagram  illustrating  an  ice 
core  (L2)  from  a  refrozen  lead  and  the 
particulate  absorption  coefficients  as  a 
function  of  wavelength  (400-700nm)  for 
the  particulates  collected  from  vertical 
sections  A  through  E  of  the  core. 


the  only  core  containing  particles  with  absorpUon  peaks  or  shoulders  in  the  430-440  nm  wavelengths  indicative  of  living 
plant  material.  Despite  these  ap(X.)'s  showing  characteristics  of  a  microalgal  component  it  is  apparent  that  particulate 
absorption  was  dominated  by  the  detritus.  The  ap(X)'s  from  an  equivelent  volume  of  the  water  column  were  2  to  100  times 
lower  than  the  apfA-l's  in  the  ice  (data  not  shown),  thus  suggesting  that  particles  are  more  concentrated  in  sea  ice  than  in  the 
water  column. 

3.5  Particles  contribution  to  K^iice  QA 


Using  the  attenuation  coefficients  of  multiyear  ice  (figure  4)  and  of  blue  first-year  ice  1.22  as  estimates  of  the 
attenuauon  coefficients  of  the  multiyear  ice  and  the  lead  ice,  respectively,  we  estimated  the  ratio  of  the  particulate 
attenuaaon  coefficient  to  the  total  sea  ice  attenuation  coefficient  (figure  10).  This  ratio  indicates  that  particulate  material 
had  their  largest  contribution  at  wavelenghts  <500  nm  where  sea  ice  generally  has  low  attenuation  coefficients!  and 
absorption  by  the  particulate  material  is  at  its  highest. 
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Figure  10.  Estimates  of  the  fraction  of  the  sea  ice  attenuation  coefficients  due  to  particulate  absorption  (F)  as  a 
function  of  wavelength  at  four  core  sampling  sites.  Upper  solid  line  =  L2,  lower  solid  line  =  LI,  open  diamonds  = 
HH2.  closed  diamonds  =  HHl.  F  calculated  according  to  equation  7.  The  attenuation  coefficients  for  HHl  and  HH2 
are  shown  in  figure  4  while  attenuation  coefficients  of  LI  and  L2  were  estimated  by  adding  the  attenuation  coefficient 
of  blue  first-year  ice21  to  the  calculated  <K<jp(A.>  for  LI  an  L2. 


3.6  EMR  absorption 


The  application  of  the  energy  absorption  model  indicates  that  particulate  matter  present  at  CEAREX  "O"  camp 
(i.e.  L2)  was  the  major  energy  absorbing  constituent  within  the  sea  ice.  The  fraction  of  the  energy  being  absorbed  by  the 
particles  within  the  ice  becomes  more  pronounced  as  snow  depth  increases  (rows  1-4,  table  1).  This  is  caused  by  the  snow 
shifting  the  transmitted  inadiance  towards  wavelengths  which  are  readily  aborbed  by  particulate  matter.  The  magnitude  of 

energy  absorption  by  both  the  sea  ice  and  the  particles  is  strongly  reliant  on  snow  depth  which  has  been  demonstrated 
previously.  1.2.3 


Figure  11.  Downwelling  spectral  irradiance  at  the  surface  of  the  ice  (measured  on  JD  99,  08(X)  hours),  and  predicted 
under  ice  downwelling  spectral  irradiances  with  and  without  300  mg  chla  m-2.  Ice  thickness  =  1.2  meters.  Snow  depth 
=  0.001  meters. 


Table  I.  Estimate  of  the  partitioning  of  energy  absorbed  by  sea  ice.  sea  ice  particulates  and  energy  passing  into  the 
water  column. 
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depth  of  ice  layer  1  and  2  respectively,  E^= 

energy  absorbed  by 

constituent  x,  where  p  =  particulate  matter,  and  i  =  ice.  The  column  headed  by  "H2O"  shows  the 
amount  of  energy  leaving  the  ice  sheet  and  passing  into  the  water  column.  %  Eabs  by  particles  is  the 
percentage  of  the  energy  absOTbed  by  the  ice  sheet  which  is  absorbed  by  the  particulate  matter. 
Incident  downwelling  irradiance  was  set  as  that  measured  at  0800  hour  on  JD  99  (figure  2  and  11) 
(152  W  m-2). 

•  ap(X)  set  as  the  <ap<X.)>  of  the  lead  ice  core  L2  (figure  9). 

-  ap(X)  set  to  zero  at  all  wavelengths. 


**  ap(X)  set  as  the  ap(A.)  of  the  particulate  matter  collected  from  section  C  of  the  HH2  ice  core  (figure  7). 

@  ap(X)'s  modeled  as  the  chla  specific  absorption  coefficient  (m2  mg  chla-^  of  a  congelation  ice  algal 
community^  multiplied  by  the  chlorophyll  concentration  in  order  to  approxunate  the  effect  of  a  rich  algal  layer 
(300  mg  m-2)  withm  the  ice.29  Figure  1 1  illustrates  the  predicted  under  ice  uradiances  with  300  mg  chla  m-2. 


4.  DISCUSSION 

It  has  been  observed  on  several  occasions  that  sediment-laden  ice  or  algal  assemblages  within  sea  ice  may  induce 
premature  ice  melting  and  alter  the  physical  structure  of  the  sea  ice  and  snow.  Cryoconite  holes,  for  example,  are  common 
in  ice  which  has  sediment  deposits  on  its  surface25  and  ablation  of  bottom  ice  occurs  earliest  where  algal  assemblaaes 
occur.26  These  phenomena  are  thought  to  be  mediated  through  the  selecuve  absorption  of  EMR  by  the  particulate  mauer 
and  re-radiation  of  the  absorbed  energy  as  heat  thus  causing  localized  melting. 

A  single  stream  radiative  transfer  model  applied  to  multiple  layers  indicates  that  particles  were  the  major 
constituents  absorbing  EMR  (400-700  nm)  within  the  interior  portions  of  sea  ice.  When  detritus  and/or  minerals  are  more 
abundant  than  observed  at  "O"  camp  or  when  miaoalgal  communities  develop  in  the  sea  ice,  they  are  predicted  to  have 
even  more  dramatic  effects  on  the  transfer  of  EMR  into  the  sea  ice  and  upper  layers  of  the  Arctic  ocean.  The  prescence  of  a 
rich  algal  layer  reduces  and  shifts  the  transmitted  spectral  irradiance  from  500  nm  to  >550nm  wavelengths  (figure  II)  (also 
see  references  5,7,  and  8).  Furthermore,  an  ice  sheet  containing  a  rich  algal  layer  in  the  bottom  portions  of  the  ice  may 
increase  the  energy  absorbed  within  the  sea  ice  by  300  to  360%  (table  1,  rows  10  and  1 1). 

Although  the  ice  cores  collected  at  "O"  camp  contained  areas  of  visually  apparent  particulate  matter,  the  majority 
of  the  ice  volume  appeared  particle  free.  Yet,  when  the  samples  were  analyzed  microscopically  and  optically  it  was  revealed 
that  all  layers  along  the  sea  iceprofile  contained  particulate  matter  which  contributed  to  the  inherent  optical  properties  of 
the  ice  sheet.  Multiple  studies25,28,28  indicate  particle-laden  sea  ice  in  the  Arctic  basin  is  common  and  is  a  significant 
carrier  for  cross-shelf  transport  of  sediments  and  continental  shelf  erosion.28  The  significance  of  particle-laden  ice 
throughout  the  Arctic  basin  is  that  biogenic  as  well  as  nonbiogenic  particulate  matter  may  play  a  significant  role  in  local 
radiative  transfer  and  therefore  biological  (e.g.  primary  production)  as  well  as  physical  processes  including  heat  and  energy 
transfer.  Their  effect  will  be  most  pronounced  in  areas  where  there  is  little  snow,  and  turbid  ice  occurs  due  to  sediment 
entrapment  and/or  miCToalgal  community  development  within  the  ice  matrix. 

Since  particulate  material  in  sea  ice  is  variable  in  both  space  and  time,  we  can  not  yet  determine  quantitatively  the 
role  particulates  have  in  determining  the  optical,  thermal,  and  biological  properties  of  Arctic  sea  ice.  Further  analysis  of  the 
factors  controlling  the  spatial  (both  vertical  and  horizontal)  and  temporal  distribution  of  terrigenous  matter,  aeolean 
particulates,  detritus,  and  igae  in  sea  ice  seems  warranted. 
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